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Abstract 

We introduce nonparaxial spatially accelerating waves whose two-dimensional transverse profiles 
propagate along semicircular trajectories while approximately preserving their shape. We derive 
these waves by considering imaginary displacements on spherical fields, leading to simple closed- 
form expressions. The structure of these waves also allows the closed- form description of pulses. 



1 



The so-called "accelerating waves" have received considerable attention due to their re- 
markable properties: they preserve their intensity profile under propagation, but the profile's 
features follow a curved path. This behavior seems at odds with Ehrenfest's theorem, which 



ree paraxial beam follows 
3[ , although more general 



in this context states that the transverse intensity centroid of any 
a straight path. The best-known example is that of Airy beams [1 
families of accelerating paraxial beams exist {4, 5|. 

The strange behavior of these beams can be understood through two observations. First, 
their transverse intensity is not integrable, so they require infinite power and their transverse 
centroid is not well defined (i.e., Ehrenfest's theorem is not violated). Second, they are 
associated with rays forming parabolic caustics, so the transverse intensity maximum is due 
to different ray bundles at different propagation distances. That is, one should not think 
of this maximum as a "physical entity" following a curved path. In real optical systems, 
the plane-wave spectrum is limited, leading to finite-power approximations to accelerating 
beams whose intensity profile is roughly preserved over some distance. While the intensity 
maximum of these beams still follows a curved path within this range, the now well-defined 
intensity centroid describes a straight line. 

Nonparaxial accelerating waves in two dimensions, whose intensity profiles are roughly 
preserved within a range of propagation distances, and whose maxima follow semicircular 
paths, were recently proposed p, lZ[. These waves result from truncating the backward- 
propagating components of fields with cylindrical symmetry. Monochromatic and pulsed 
waves following other nonparaxial non-circular paths have also been implemented jsl,]^]. 

In this Letter, we extend these ideas to three-dimensional waves by starting with fields 
with rotational symmetry. Our solutions have a two-dimensional transverse intensity profile 
that is approximately invariant under propagation and whose maxima follow circular paths. 



To construct these waves we apply the concept of imaginary displacements 10] . By this 



means, unlike the solutions generated by imposing an abrupt exit pupil js, |7|, we are able 
to give simple close-form solutions. 

We begin by considering two-dimensional waves, i.e., solutions to the two-dimensional 



Helmholtz equation. One such solution is the Bessel field 



Hi 



kVx 2 + 



exp [ira arctan(2;, x)] 



where u is a constant, m is an integer (assumed here to be positive), J m is a Bessel function 
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of the first kind, k is the wavenumber, (x, z) are Cartesian coordinates where z is regarded 
as the propagation direction, and arctan(cos 0, sin <ft) = <p. This field has an intensity profile 
that is exactly preserved not over lines of constant z but over lines crossing the origin, and 
its maximum follows a circular caustic of radius m/k. A section of this profile is shown in 
Fig. QJa). For large m, the radial profile of this field approaches the profile of a paraxial 
Airy beam over a finite radial region near m/k. In fact, accurate approximations to these 



fields can be written in terms of Airy functions even for modest values of m 
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FIG. 1. (Color online) Intensities over a 90 x 80 rectangle (in units of k^ 1 ) of m = 40 Bessel fields: 
(a) unapertured, apertured with (b) a = tt and (c) a = n/2, and apodized with (d) q = 1.88 and 
(e) q = 3. The insets indicate the power angular spectrum. 



Notice, though, that Bessel fields are composed of plane waves traveling in all directions: 

/' TT 

ip m (x,z) = I A m ((j)) exp[ifc(xsin0 + zcos0)] d(/>, (2) 



where <fi is the plane waves' propagation angle from the z axis, and A m ((p) = uq exp(im</>)/27r 
is the angular spectrum. In practice, it is difficult to generate fields whose plane-wave 
spectrum extends over all directions, so fields with spectra confined to forward propagation, 
i.e., A ((f)) = for cos<ft < 1, are usually considered. For this purpose, Kaminer et al. 
6| proposed limiting the integration in Eq. ([2]) to the forward semicircle G [— 7r/2, 7r/2], 
which results in a forward-propagating wave with accelerating characteristics. Smaller ranges 
4> G [—a/2, a/2] for < a < n can also be used (and would be easier to implement in 



practice). The intensities of these fields for a = tt and tt/2, are shown in Figs. QJb), EJc), 
respectively. Notice that the intensity maximum propagates approximately along a circular 
path until its angle with respect to z is slightly below a/2. However, for any a 7^ 2ir, the 
fields must be evaluated numerically, and as shown in Fig. mb)-QJc), the hard integration 
limits introduce appreciable rippling in the intensity. 

To avoid these issues, we propose to use instead a Gaussian-like apodization factor 

k 



equivalent to multiplying by 
10] . The resulting field is given 



exp[/cg(cos0 — 1)] in the angular spectrum. This is forma 
exp(— q) and performing an imaginary shift z — > z — iq/k 
simply by 

i/) m {x, z; q) = exp(-g) ip m (x,z - iq/k). (3) 

The intensities of these fields are shown in Figs. QJd)-flJe) for q = 1.88 and 3. Note that 
the intensity maximum decreases smoothly away from z — 0, and follows a circular path 
until this path forms an angle of about 2 &rccos[Ii(q) / Io(q)} with respect to z, with I n being 
the modified Bessel function of the first kind. While these fields contain some counter- 
propagating components, these represent only a fraction [1 — Lo(2q) / Io(2q)}/2 of the total 
power, where L n is the modified Struve function of the first kind. For the cases in Figs. H(d)- 
[He), this power fraction is about 1% and 0.082%, respectively. 

We extend these ideas to create three-dimensional accelerating waves by using multipoles 
(separable in spherical coordinates), whose plane- wave spectra are spherical harmonics Y^ m , 
aligned with the y axis for convenience: 

A™(r) = J >7 m (#,0)exp(i&r-u) dfi 
= 4iri l ji(k^/x 2 + y 2 + z 2 ) 

x arctan (y, V x 2 + z 2 j , arctan(z, x) , (4) 

where r = (x,y,z), ji is a spherical Bessel function of the first kind, and the integral is 
over the unit vector u = (sin # sin 0, cos#,sin#cos0). We then define spherical waves with 
imaginary displacements as 

K (r; v) = u o eM-q)K+n y, * - iq/k) , (5) 

for n > 0. For q = (no apodization) these fields have exact rotational symmetry about 
the y axis. For q>2, on the other hand, the counter-propagating components are largely 



suppressed, and for m 3> n, these fields can be regarded as nonparaxial versions of Airy- 
Hermite- Gaussian beams, whose amplitude profile in y is a Hermite-Gauss function of order 
n (i.e., n specifies the number of zeroes in the y direction), while in the xz plane they 
behave as Airy beams. The intensity of some of these fields is shown in Figs. [2](a)-I21(h). Of 
course, instead of apodizing through an imaginary displacement, one can define apertured 
versions where the integral in Eq. (J4]) is limited to a forward-propagating section of the sphere 
[sin 9 cos <p > cos(a/2)] and must be evaluated numerically. The case of half-spherical waves 
(a = n) is shown in Figs. E]^i)-|2]d). Whether apodized or apertured, these spherical waves 
present an intensity maximum (or maxima for n > 0) that follows a circular path of radius 
slightly larger than m/k, cf. Figs. [2](a,e,i). 
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FIG. 2. (Color online) Intensities over sections of size 180 x 156 (in units of A;" 1 ) of the z = and 
x = planes for (a-h,m) the waves in Eq. ([5]), (i-1) truncated spherical waves with a = tt, (n) the 
waves in Eq. © , and (o-q) the pulses in Eq. ([8j) (movie online) . The cross-section of the caustic is 
shown in (m) and (n). 



It is seen from Fig. [2]that, for larger n, the field's transverse intensity profile is noticeably 
curved. This is because the caustic, whose cross-section is shown as a white line in Fig.^m), 
is composed of segments of a sphere of radius (m + n)/k and a (two-sided) cone of half-angle 
arcsin [m/ {m + n)] around the y axis, both centered at the origin. The intensity then presents 



two cusps, one at each intersection of the sphere and the cone, and therefore resembles two 
joined Airy beams. One of the cusps can be suppressed, e.g., by combining three of these 
fields as 

K - \ - K-i) . (6) 

The intensity of this field over the z = plane is shown in Fig. [2]^n). Note that it indeed 
resembles a single Airy beam whose cusp at the z = plane is approximately at y — 
a/ (2m + n)n, x = — m, but where one of the caustic sheets is noticeably curved. For n ~ 
(\/2 — l)m 3> 1, this curvature becomes less appreciable and the field increasingly resembles 
an Airy beam. 

The fields described so far are monochromatic. We now define closed-form pulsed solu- 
tions to the wave equation that follow similar paths. This is possible because in Eq. @, 
k appears only as a factor in the argument of ji, which can be written as a finite sum of 
powers times exponentials of its argument: 

. ^ {21 -p)\ [(-ra)Pexp(ra) - (ra)Pexp(-ra)] 

3l{U> X-P+Hu^plil-p)] ' {> 

Substituting u = ur/c (where ui is the frequency and c the speed of light) in Eq. flTJ) and 
multiplying by u l+l leads to an expression where each term includes u as a non-negative 
power times an exponential. By using this result together with Eq. (j4j) and applying standard 
Fourier relations, we define the following pulsed solutions 

P n m M) = J /»a; m+ " +1 A™ +n (r)exp(-ia;t)du; 

. m+n+l 

2r J 

(2m + 2n — p)\ ( 2r s 



x 



E 



^ p\ (m + n — p) ! \ c 



x [(-V P f ip) (t + l)-f {p) (t-l)}, (8) 



where f(co) is an arbitrary spectrum, f(t) is its inverse Fourier transform, is the pth. 
derivative of /, and (r, 9, <j)) are spherical coordinates with respect to the y axis. Figures^o)- 
|2^q) show cross sections at three times for such a pulse for f(t) = exp(— t 2 A 2 /2 — iuot) [for 
which can be found in closed form in terms of Hermite polynomials] with A = uq/10, 
m = 50, n = 0, and using an imaginary displacement in z of size 1.88c/u;o. As in the 
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monochromatic case, for larger values of n these pulses display two cusps, one of which can 



be eliminated throug 
unlike the pulses in 



i the combination P™ — (Pfi^i + P%L 1 )/2. It is worth noting that, 
9|, these pulses do not run along the caustic, but rather inhabit it 
simultaneously at a given range of times. 

In conclusion, we proposed nonparaxial accelerating fields given by simple closed-form 
expressions. While these fields are scalar, polarization can be easily incorporated by either 
using vector multipoles or applying suitable operators. Also, the paths followed by the 
maxima are circular due to the use of solutions separable in spherical coordinates. Using 
other coordinate systems would result in other shapes of the caustic sheets. 
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